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Abstract. Proposed by Wolfgang Pauli more than 80 years ago, the exclusion principle has been proven to have 
a far-reaching consequence, from femtoscopic world to macroscopic, super-dense, and fully relativistic physics. 
Starting from this principle, we discuss two interesting research topics, which have currently drawn considerable 
attention in the nuclear- and astrophysics communities; the hypernuclear and neutron star physics. Special attention 
is given to the electromagnetic production of the hypertriton and the consequences of the neutrino electromagnetic 
properties in dense matter. We also touch on the new arena which could also be fascinating for physicists; the 
econophysics. 

1. Pauli Exclusion Principle and Its Consequence 

1.1. Historical Background 

Most of us started to be acquainted with the Pauli exclusion principle at high-school, when we learned 
about atomic shells in the chemistry or in the physics course. Probably, many of us are ignorant that the 
discovery of the exclusion principle is in fact intimately related to the conception of the electron spin 
and the development of quantum mechanics. If we try to comprehend the motivation of Wolfgang Pauli 



behind this great discovery, for which he was awarded with the Nobel Prize in 1945, then we will find 
that it was to answer the problem of explaining the observed atomic spectra and the "anomalous Zeeman 
effect" [1]. At that time these transitions were particularly confusing, because some of them (e.g., the 
Sodium atom) exhibit doublet transitions, whereas some others (e.g., the Magnesium atom) show either 
[ — . | singlet or triplet transitions [2]. Furthermore, the anomalous Zeeman effect yields a doublet splitting, 

contrary to the what was hitherto known by the spectroscopists as triplet splitting. Another disturbing 
fact to Pauli was the question as to why all electrons for an atom in its ground state were not bound in the 
innermost shell, a problem which had been emphasized even by Niels Bohr as his fundamental problem. 

Before 1924, Arnold Sommerfeld, Alfred Lande, and Pauli himself realized that the answer to these 
problems was the existence of a new (the fourth, after n,£ and m, to be more precise) quantum number. 
Sommerfeld and Lande considered this number as an intrinsic property of the atomic core, without 
specifying whether it is the nucleus or the nucleus plus inner electrons. This point of view was then 
considered by Pauli as rather orthodox. 

In 1924 Pauli published an objection to this argument [4]. Instead of considering the properties of the 
atomic core, Pauli proposed a new quantum property of the electron, which he called "two-valuedness not 
describable classically." Using this idea Pauli was able to reduce the complicated numbers of electrons 
in closed subgroups to the simple one, once this fourth quantum property is introduced. Nevertheless, 
with the exception of few famous spectroscopists, nobody could easily follow this idea, since there was 
no analogy in terms of mechanical model. This difficulty was then remedied by the concept of electron 
spin proposed by George E. Uhlenbeck and Samuel A. Goudsmit, who tried to explain the experimental 
result of Otto Stern and Walther Gerlach [5]. It was said that after some discussions with Lorentz, they 
tried to withdraw their paper, but it was too late and the paper was published! With this concept, i.e. by 
assuming that the spin quantum number of electron as A, it is then possible to understand the doublet 
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Fig. 1. Wolfgang Pauli (1900-1958), the 1945 Nobel laureate for "the discovery of the exclusion principle" [3]. 

splitting in the anomalous Zeeman effect [6]. Since then the idea of exclusion principle was always 
closely connected with the spin. 

Pauli formulated his exclusion principle comprehensively in 1925 [7]. At about the same time, sev- 
eral milestones in quantum mechanics were achieved. The most important and relevant one to the present 
discussion is the formulation of the quantum mechanics by Werner Heisenberg [8]. According to Heisen- 
berg, quantum mechanics leads to qualitatively different conclusions for distinguishable particles than 
for indistinguishable particles. It was found that for indistinguishable particles the wave function must 
be either symmetric or antisymmetric. Pauli considered the antisymmetric one as the appropriate wave 
mechanical formulation of the exclusion principle, since an antisymmetric wave function containing two 
particles in the same state is equal to zero. 

The first indication that proton has spin ^ was obtained from the investigation of the anomaly in the 
specific heat of the molecular hydrogen in 1927 [9]. Later on, this was confirmed by Stern by using a 
Stern-Gerlach type experiment. In 1932, the neutron was discovered by Frederic and Irene Joliot-Curie 
along with James Chadwick. Two years later, by dissociating deuteron with gamma ray, y+d —>n + p, 
Chadwick and Maurice Goldhaber were able to show that a deuteron was made of a proton and a neutron. 
From its spectrum it was known that the deuteron has spin 1. Hence, we know that the neutron has also 
spin ^ and both proton and neutron obey the Pauli exclusion principle. 

Although for us the exclusion principle is quite clear and understandable, Pauli himself was dis- 
appointed with his finding, since he said that he was unable to give a logical reason or to deduce the 
principle from a more general assumption. He was hoping that the development of quantum mechanics 
will also rigorously deduce the principle, but he eventually found that it was not the case [1]. 



1.2. The Degeneracy Force 

Let us begin with a simple gedanken experiment, i.e., we drop a number of fermions into a square well 
potential. We repeat the same experiment for a number of bosons. As we expected, the result is shown 
in Fig. [2] However, to understand why the fermions sitting above the ground state do not collapse to the 
lowest state as in the case of bosons, we need to understand the degeneracy force, a kind of force which 
maintains such formation. 

Consider two particles a and b located at x\ and The wave function of such a system can be 
written as [10] 

x P(x 1 ,X 2 ) = \f a (x\)\fh(X2) , (1) 

where both \\r a and \|f& are normalized to unity. Let us now calculate the expectation value of the "dis- 
tance" between the two particles, i.e., 

<(jci ~X2?) = (4) + (4) -2<jcijc 2 ) • (2) 
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Fig. 2. If nine fermions and nine bosons are "dropped" into square well potentials, they will organize themselves 
in quite different ways. This happens because the fermions always obey the Pauli exclusion principle. 



Using the wave function given in Eq. CD we can calculate the first term on the r.h.s. of Eq. (0, 

{x\) = j dxix\\y a {xi)\ 2 J dx 2 \yb{x2)\ 2 = J dx x 2 \y a (x)\ 2 = (x 2 ) a , (3) 

N v ' 

= i 

where the subscript a indicates that the expectation value is calculated using the wave function \|/ a (x). 
By following the same procedure we can easily prove that {x 2 ) = (x 2 )b and {x\X 2 ) = (x) a (x)b- Therefore, 
the "distance " between two distinguishable particles is found to be 

-x 2 ) 2 )dis = (x 2 ) a + (x 2 ) b -2(x) a (x) h . (4) 

On the other hand, the wave function for two indistinguishable boson or fermion is given by 

x ¥±(x U X 2 ) = -^={\\f a (x l )\\fb(x 2 )±\\f a (x 2 )\\f h (x l )} , (5) 

where the + and — signs indicate the symmetric (boson) and anti-symmetric (fermion) systems. Using 
Eq. ((D we can recalculate Eq. (0]) and obtain 

((x l -x 2 ) 2 )± = (x 2 ) a + (x 2 )b-2(x) a (x)bT2\(x) ah \ 2 , (6) 

where the last term, (x) a b = J dx^*(x)x^b(x), clearly depends upon the overlap between the wave 
functions y a (x) and 

By comparing Eqs. (@]) and (© we can conclude that 

{{xi-x 2 ) 2 )± = ((xi-x 2 ) 2 ) dis ^2\{x) ah \ 2 . (7) 

This equation is extremely important. Although derived in a quite simple fashion, it shows that two 
fermions tend to be farther apart than two bosons. This is known as the exchange force (or degeneracy 
force), which produces a "force of attraction" between identical bosons and a "force of repulsion" be- 
tween identical fermions. Note that this is not really the force we are familiar with, since there is no 
mediator which pull or push the particles, instead it is solely a consequence of the (anti)symmetrization 
required by quantum mechanics. 

After we know that identical fermions exhibit "repulsive force", our next question is certainly "how 
repulsive is this force?" To answer this question we use the electron gas model, since in this model 
electrons do not interact via Coulomb force. Furthermore, the electron gas model is a simple model 
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Fig. 3. A A-hyperon can freely move in the nuclear energy states, without being blocked by the Pauli principle. 



often used in many standard quantum mechanics books [24]. In this model, for N number of electrons, 
each with a mass of m, occupying a volume V, the corresponding total energy can be written as 



h 2 % 3 /3N\ 5/3 



WmV 2 / 3 \n J 



- 



(8) 



Using the density p = N/V we can calculate the "pressure" produced by the degeneracy force of N 
electrons and find that 

dE tnt _ fcV/3p\ 5/3 



Pd ^ = "air = T5^ UJ • (9) 

To get a better feeling on this pressure, let us consider two "chargeless" electrons separated with a 
distance of 1 A. The particle density in this case is p 10 30 m~ 3 . Substituting this value into Eq. (© we 
obtain the degeneracy pressure is pdeg. ~ 10 14 N/m 2 . One may suspect that such a huge pressure would 
be meaningless in an atomic scale, where the typical "area" is of the order of 10~ 20 m 2 and therefore the 
corresponding force is of the order of 10~ 6 N. Nevertheless, one should not forget that at this scale the 
Coulomb force between two electrons is also of the order of 10~ 8 N. The problem, however, turns out to 
be dramatically interesting on the macroscopic scale where we do business with an area greater than 1 
km 2 or 10 6 m 2 . We will encounter this interesting problem when we discuss the neutron star, but first we 
will introduce the physics of hypernucleus. 



2. Introducing a A Hyperon in the Nucleus; Hypernuclear Physics 

2.1. Pauli Principle in the Nucleus 

Let us now imagine that we can implant a A hyperon inside a normal nucleus. Since the A hyperon is 
distinguishable from the nucleons (the protons and neutrons which build up the nucleus), it will certainly 
be free from Pauli blocking (see Fig. 0. Such a nucleus is called a hypernucleus. Hypernucleus is 
interesting because many new information about the nuclear structure can be obtained by studying the 
way the A hyperon localizes inside the nucleus. Another interesting aspect comes from the fact that a 
free A hyperon is unstable, it decays in about 10~ 10 s primarily to 

J p + %~ (branching fraction : 63.9%) 
~* \ n + %° (branching fraction : 35.8%) ' ( ' 

via weak interaction. However, if the A is embedded in a nucleus, since the low energy nucleon states in 
the nucleus are already filled, Pauli principle will prevent it from decaying in the usual way as in Eq. (flOl ). 
Therefore, there must be another mechanism which eventually allows the hypernucleus to decay via a 
weak process. 




Fig. 4. The first observed hypernucleus. A highly energetic proton p strikes a nucleus at point A in an emulsion, 
which immediately breaks up into a number of nuclear fragments. All fragments are eventually stopped in the 
emulsion except the one marked with f. This fragment decays at point B into three charged particles and a neutron. 
The neutron is not seen in this figure, since it is chargeless. Original figure was taken from Ref. [12]. 

The first hypernucleus was observed by Marian Danysz and Jerzy Pniewski, two Polish physicists, 
in 1953. They used photographic emulsions exposed to the cosmic radiation. These emulsions recorded 
a break-up event illustrated in Fig. [4] We can understand this event as follows. A highly energetic proton 
coming from cosmic rays hits a heavy nucleus in the emulsion, causing it to break up into a number of 
fragments. One of them (marked with f) decays into three charged particles at point B. This decay cannot 
be interpreted as a normal decay of a highly excited nucleus, since from the track shown in the figure 
the three particles should have lived for about 10~ 12 s, whereas the normal decay of an excited nucleus 
would take place in about 10~ 22 s. Presumably, the event can be interpreted as the decay of the ^Be 
hypernucleus into three charged particles [11], 

^Be -> 3 He + 4 He + n + 7t~ . (11) 

The first hypernucleus was studied in a photographic emulsion exposed to cosmic rays. Nowadays, 
hypernuclei can be produced in laboratories by using, e.g., hadronic reactions 

K + A Z -» A A Z + 71" , 7C+ + A Z -» \Z + K + , (12) 

or through electromagnetic interactions 

y + A Z -» \Z + K + , e + A Z -» e' + A K Z + K+ . (13) 

The advantage of using the first reactions, Eq. (fl2l . is that the cross section is large and therefore the cor- 
responding experiments are relatively easy to perform. However, the electromagnetic reactions, Eq. (fT3T >. 
are also invaluable for the study of A hypernuclei. Since a large momentum is transferred to the recoiled 
hypernucleus, it will then populate high spin states and, therefore, is very suitable for investigating deeply 
bound states of a A hyperon. Furthermore, both spin-flip and non-spin-flip amplitudes are excited with 
significant cross sections, in contrast to the hadronic reactions. In addition, since a proton is converted 
to a A hyperon, hypernuclei that are not accessible by hadronic reactions can be produced. Finally, a 
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Fig. 5. (Left) Photoproduction of a hypertriton. (Right) Photoproduction of a A hyperon which is called elemen- 
tary reaction. The amplitude obtained from Feynman diagrams for this reaction is called elementary operator. 

comparison between hadronic and electromagnetic reactions allows us to study charge symmetry in the 
Ahypernuclei [13]. 

The simplest example is photoproduction of a hypertriton, the lightest hypernucleus. A hypertriton 
consists of two nucleons and one A hyperon. Thus, it can be produced by bombarding a 3 He target with 
(real or virtual) photon, as shown in Fig. [5] From this figure it is clear that one of the proton inside 
the nucleus is converted into a A hyperon along with a kaon through an elementary process. Therefore, 
to theoretically calculate the hypertriton photoproduction one needs information about the elementary 
reaction. 

2.2. Elementary Reaction 

The elementary operator can be modeled in terms of Feynman diagrams at tree-level as shown in Fig. [6] 
Although higher-order corrections are neglected, most of the hadronic and electromagnetic couplings 
in these diagrams are in fact not well known. Fortunately, experimental data for the elementary process 
y+ p — > K + + A have been available since more than 40 years ago. Thus, all unknown coupling constants 
can be extracted by fitting the calculated observables to these data. 

Very recently a large number of data have been published by the SAPHIR [14] and CLAS [15] 
collaborations. These data cover the energy range from reaction threshold up to W = 2.5 GeV, angular 
distributions from forward to backward kaon angles, and are significantly more accurate than the previous 
ones [16]. A comparison between theoretical calculations with small part of these data is shown in Fig. [7] 
To calculate the differential cross section we have used 

do \q[ m p E A 2 

d£l \k\ 32tu 2 W 2 V ' ' 

where q and k are the kaon and photon momenta, respectively, whereas W is the total cm. energy. The 
matrix element is obtained from the Feynman diagrams shown in Fig. [6] and is summed over the 
possible photon polarization £. 

A quick glance to Fig.|7]can immediately convince us that the new data have a problem of the mutual 
consistency. To study the physics consequence of this problem, Refs. [17] and [18] have performed 
different fits to several data sets. In Fig. |7] we show only two of them, i.e., fit to the new SAPHIR data 
(Fit 1) and fit to the new CLAS data (Fit 2). The new LEPS data are found to be more consistent with 
the CLAS ones. 

Reference [18] has pointed out that the use of SAPHIR and CLAS data, individually or simultane- 
ously, leads to quite different resonance parameters. This could lead to different conclusions on "missing 
resonances", the resonances which have been predicted by constituent quark models but are intangible 
\o% + N ^K + N reactions that are used by the Particle Data Group (PDG) to extract their properties. 
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Fig. 6. Elementary operator of KA photoproduction on the nucleon can be modeled by several Feynman diagrams 
at tree-level. Diagrams (a), (b), and (c) correspond to the s-, u-, and f -channel, respectively, whereas diagram (d) 
represents the contact term which is required to maintain gauge invariance of the process. 



In view of this, the elementary process of kaon production can also be used as a complimentary tool 
for investigation of the "missing resonances". As an example, it is worth mentioning the case of "miss- 
ing" Di3(1895) which was observed by Ref. [19] in the y+ p — ► K + + A channel. It was found after 
the previous SAPHIR data [16] showed a clear peak around W = 1900 MeV in the differential and total 
cross sections. Subsequently, it was shown in Ref. [20] that the peak could also be equally well repro- 
duced by including a Pi3(1950) resonance. However, most of analyses based on the isobar model after 
that confirmed that including the Di3(1895) significantly improves the agreement between theoretical 
prediction and experimental data. Recent partial-wave and coupled-channels analyses confirmed that the 
peak originates from a D13 resonance with a mass around 1900 MeV [21]. Indeed, the same conclusion 
has also been drawn from the recent multipole study [19], that the peak at W ~ 1900 MeV originates 
from a D13 resonance with the extracted mass equals 1936 MeV if SAPHIR data were used or 1915 MeV 
if CLAS data were used. Thus, future measurements such as the one planned at MAMI in Mainz are, 
therefore, expected to remedy this unfortunate situation. 



2.3. Electromagnetic Production of the Hypertriton 

Theoretical calculation of photoproduction of the hypertriton can be done by making use of the left 
diagram of Fig. [5] The formula for calculating the hypernuclear cross section is similar to the one for 
calculating the elementary cross sections [cf. Eq. (fT4V l. i.e., 

da jgj M3 He £3 H 

da |*| L m L 2 1 fil > ^ 

except now the nuclear amplitude squared is averaged over the initial and final nuclear spins. The 
nuclear amplitude is somewhat complicated in this case. It consists of the elementary part, integrated 
over the possible initial and final nucleon momenta and weighted by the nuclear wave functions, and 
some angular parts in terms of 6 j and 9j coefficients. Explicitly, it is given by 



* ^ an 1 11 A 111 * 
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Fig. 7. (Left) Comparison between experimental total cross sections with the predictions of Fit 1 (fit to the SAPHIR 
data) and Fit 2 (fit to the CLAS data). (Right) Differential cross sections obtained from Fit 1 and Fit 2 as a function 
of the total cm. energy for different kaon angles. Solid circles, squares, and triangles represent experimental data 
from the CLAS, SAPHIR, and LEPS collaborations, respectively [18]. 



J d 3 qp 2 dp <p A (q f ) Vf\p) UP,l) [YW(f )®K(" 



(A) 



where the amplitude is related to the elementary amplitude given in Eq. (fT4l ) by 



L + ia-K = £ (-ifn [oW®kW 



«=o,i 



(0) 



(16) 



(17) 



with a^ 1 ' is the Pauli matrices. 

The result is shown in Fig. [U where we compare the predicted cross sections of hypertriton pro- 
duction with that of elementary reaction. For the hypernuclear production we have performed an exact 
calculation using Eq. (fT5l ) and an approximated one by using the so-called "factorization" approximation. 
In the latter, the cross section is calculated by means of 



(18) 



elementary 



where Wa is a kinematical factor and the "form factor" F(Q) is given by the overlap integral of the two 
wave functions, 



HQ) 



(19) 



Note that in Eq. ( fl9l ), p and q denote the momenta of the interacting nucleon and the spectators inside 
the 3 He, respectively, whereas Q is the momentum transfer (see Fig. [5]). 

The nuclear cross section at forward angles is smaller than that of elementary kaon production by two 
orders of magnitude. As 6^ m increases, the cross section drops quickly, since the nuclear momentum 
transfer increases as a function of 6^ m . Since this cross section is in fact very small, photoproduction of 
hypertriton provides a great challenge to experimentalist. Figure [8] also shows the significant difference 
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Fig. 8. Differential cross section for kaon photoproduction off the proton (dotted line) and 3 He (dashed and solid 
lines) as a function of kaon angle. The dashed line shows the "factorization" approximation for production off 
the 3 He, whereas the solid line represents the exact calculation by using S- waves [22]. Solid circle represents the 
CLAS datum at this kinematic [15]. 



between the cross sections calculated with the approximation of Eq. (TT8T ) and the full result obtained from 
Eq. ( fT5T >. This discrepancy is due to the "factorization" approximation, since in the exact calculation 
both spin-independent and spin-dependent amplitudes are integrated over the internal momentum and 
weighted by the two wave functions. This also explains why more structures appear in the cross section 
of the exact calculation. 

Very recently, Dohrmann et al. [23] measured the electroproduction process on both 3 He and 4 He 
targets. Although their experiment yielded only three data points in each target, this provides the first 
evidence on the possibility to produce hypernuclei via electromagnetic process. Our investigation of this 
process is still in progress. To compare our photoproduction results with these new data, we multiply 
Eq. ([TBI with the ratio between the electro- and photoproduction cross sections of the elementary process. 
The results for both targets are displayed in Fig. |9l where we see that a fair agreement is obtained in the 
case of the hypertriton production and a significant discrepancy appears for the 4 He target. This clearly 
shows that our simple approach urgently requires substantial improvement. 
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Fig. 10. A solid sphere with radius R and density p. The mass of the shell with the thickness dr is Anpndr, 
whereas the mass of the sphere enclosed by this shell is 47tpr 3 /3. 



3. Moving to Neutron Star 

3.1. Simple Electron and Neutron Stars 

Let us assume that a star has burned all of its "fuel" and the gravitational force starts to contract the star 
until nothing but the electron degeneracy pressure tries to save the star from a gravitational collapse. To 
calculate the gravitational pressure, let us consider a solid sphere depicted in Fig.QIjl From this figure, it 
is straightforward to derive the gravitational potential of the shell lying between r and dr, which in this 
case can be written ad3 

_ _G(^)(4*p^ f ) _ _i^ r 4 
r z 3 

where p is the density of the sphere and G is the universal gravitational constant. By integrating over the 
entire radius of the sphere we obtain the total potential energy, 

16TC 2 p 2 G f R 4 167i 2 p 2 /? 5 G 3 2 /47i\ 1/3 1/3 

where M is the total mass of the sphere and V is its volume. The gravitational pressure is obtained by 
differentiating Eq. (f2TT > to the volume V, 

Pgrav. = = -\{^y 3 G(Nm n fV-^, (22) 

where we have assumed that the sphere contains N neutron each with a mass of m n , i.e., M = Nm n . 
The degeneracy pressure is obtained from Eq. ©, which can be written as 



Pdeg. 



15m P V 71 



^ V- 5 / 3 , (23) 



with N e and ra e are the number and the mass of electron. 

If there is an equilibrium between the two forces, then we have 

1 /4tt\ 1/3 , n2 4/ , hV /3N e \ 5/3 

«(WV^ = T^(^) y- 5 '^ (24) 

where and N e are the number of neutron and electron in the star, m„ and m e are the neutron and electron 
masses, respectively. Solving for the star radius, we obtain 

R - f 8l7t2 V /B A 2 " 3 jy-i/3 (25) 
\ 128 / Gm e ml 



'For the quantum-mechanical derivation of the neutron star radius, we follow Gasiorowicz's book [24]. 
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For a star of one solar mass, we have N « 1.2 X 10 57 and assuming N e = N/2 we obtain 

R e , ps 1.1 x 10 4 km . (26) 

Such a star, for which its collapse is halted by electron degeneracy, is called white dwarf. 

However, if we started with a heavier star then the gravitational force would be too strong and the 
electron degeneracy force could not overcome it. As a consequence, electrons and protons would start 
the process 

e~ + p -» v e + n . (27) 

This process indicates that all atoms in the star start to collapse. Since the neutrino interacts very seldom 
with other particles, it will exit the star and we are left with a pure neutron star. To calculate the radius 
of this star we can use Eq. (l25l) and substitute the mass of electron with the mass of neutron, i.e., 

If we started with a star with mass equals to the solar massJl we will then end up with a neutron star with 
R ns 10 km. This result should be compared with the solar radius which is approximately 7 x 10 8 km. 

3.2. A Slightly More Realistic Neutron Star 

The neutron star described in the previous subsection is clearly far from realistic. First of all, because a 
real neutron star contains not only neutron. Why? Let us consider again the process shown in Eq. (|27T ). 
Given that a free neutron will decay in about 15 minutes via 

n -> p + e~ + v e , (29) 

certain amount of the protons should be present in the neutron star to prevent this weak decay. Most of 
the energies in the reaction of Eq. d29l) is carried away by the anti-neutrino and, therefore, the protons 
can only occupy the low-energy levels. Once all these levels are filled, Pauli exclusion principle prevents 
the decay from taking place. On the other hand, the appearance of these protons must be accompanied 
by the presence of the same amount of electrons, in order to make the charge of the star neutral. 

Inside the neutron star the pressure is increasing. As a consequence, other particles such as /u, n, and 
K mesons may appear. Going deeper and deeper we will meet exotic particles, such as A,£,E, and A, or 
even the up, down, and strange quarks. This can be nicely seen in Fig.QT] 

Besides that, since the pressure is extremely strong, relativistic corrections start to play a significant 
role. To understand this correction, let us go back to the solid sphere shown in Fig.[l0] The gravitational 
force between the mass of the shell (dm) and the mass enclosed by it (m) can be written as 

Gmdm Gm(r)p(r)A(r)dr 
at = -J — — -j ' ( -- 5U - ) 

where p(r) is the density of the shell and A(r) is the area covered by the shell. By using dp = dF /A, we 
obtain 

dp_ _ Gp(r)m(r) 
dr r 2 



2 Only for example. This is given in order to be in line with the discussion of Gasiorowicz. A more consistent calculation 
indicates that this could be true only if we started with a star with the mass greater than 1.4 of the solar mass. 
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R ~ 10 km 



Fig. 11. The structure and phase of subatomic matter inside the neutron star as predicted by theory. Also shown in 
this figure, the density of the crust layers [28] (Courtesy of Prof. Fridolin Weber, San Diego State University). 



We also note that the mass m(r) in Eq. (f3TT > is related to the density p(r) by 



dm 
dr 



4nr 2 p(r) . 



(32) 



Equations (f3TT > and d32l provide the tools to investigate a Newtonian star, i.e., a star for which the mass 
is not so large and therefore no relativistic correction is needed. These equations can be integrated self 
consistently, from r = (where p = po and p = p c ) up to r = R (where p = 0), to obtain the pressure 
p(r) and, therefore, the radius R and the mass M of the star. To perform this integration, one needs the 
information of the matter density p(r) in terms of the pressure p(r). This relationship is well known as 
the equation of state (EOS) of the matter building up the star. 

There is a certain class of stars for which the corresponding EOS is quite simple and takes the form 



(33) 



where K is the proportionality constant which depends on the entropy per nucleon and chemical compo- 
sition. Any star fulfilling Eq. (l33l) is called a polytrope. The case of polytrope is of interest since there 
are analytic solutions in terms of Lane-Emden function 8(^) [25]. In this case the radius of the star is 
given by 



R 



Ky 



I 1/2 



4jcG(y-1) 



Pc Si 



(34) 



whereas its mass is given by 



M 



47ip 



(3y-4)/2 



Ky 



4nG(y-\) 



3/2 



3 



(35) 



where both £i and 0'(^i) are functions of the polytropic index yand tabulated in Ref. [25]. By eliminat- 
ing p c in Eqs. (l34l and (I35T ). it can be shown that the relation between mass and radius of this type of star 
is independent of the choice of p c . 
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r ( km ) 

Fig. 12. The pressure and mass of a neutron star calculated from Eqs. (f36t and d32i > as a function of the distance 
from its center. In obtaining this result we have used p c = 5po, where po = 0.17 fm~ 3 , the normal nuclear density 
[29]. 



Most of the known stars can be sufficiently described as Newtonian ones. Besides that, the compact 
stars known as white dwarfs also fall into this category. 

For more massive stars, such as neutron stars, the relativistic correction is needed. This can be 
performed by multiplying Eq. (I3TI ) with three correction factors, 

dp Gp(r)m(r) 
dr r 2 

relativistic corrections 

while maintaining the differential equation for m(f) [Eq. (l32l)1. Equation (l36l) is known as the Tolman- 
Oppenheimer-Volkoff (TOV) equation [26]. In this equation the first two factors in square brackets are 
due to special relativity corrections of order v 2 /c 2 . It is obvious that these factors go to zero in the limit of 
c — > oo. The third factor is a general relativity correction. It is also clear that this factor can be neglected 
in the non-relativistic case. It is obvious from Eq. (l36l > that an analytic solution of the TOV equation is 
far from possible. The discussion of the solution to this equation by using an advanced EOS is clearly 
beyond the scope of this paper. The interested reader can consult Ref. [27] (as well as references therein) 
for the latest review on this topic. As an example of the numerical solution to the coupled equations (l32l 
and (1361 we show the pressure and the mass of neutron star as a function of the distance from the center of 
the star in Fig. [12] Here, a simple form of the energy density of nuclear matter consisting a compressional 
term, a symmetry term, and binding energy term, has been used for the corresponding EOS [29]. 

3.3. Neutrino Mean Free Path in Neutron Stars 

In discussing the formation of a neutron star, i.e., Eqs. (l27l) and d29l ), we understand that both neutrino and 
anti-neutrino immediately leave the neutron star as soon as they have been created. Moreover, it is also 
comprehended that most of the energies coming from the two reactions are carried away by neutrino and 
anti-neutrino. This has an important consequence in the cooling process of a neutron star, i.e., whereas a 
conventional body dissipates its energy through photon emissions (electromagnetic radiations), a neutron 
star cools down via emitting neutrinos and anti-neutrinos. 

Therefore, understanding of the interaction between neutrino (as well as anti-neutrino, from now on 
we will only mention neutrino for the sake of brevity) and neutron star is crucial in the investigation of the 
neutron star formation. Since the neutron star density is extremely high (see Fig.fTTb. such interaction is 
more probable than in the conventional matter. In fact, the neutrino mean free path (the distance between 
two interactions, from now on called NMFP) could in certain cases be smaller than the diameter of the 
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P( r ) 
p(r)c 2 
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471 r p(r) 
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Fig. 13. (Left) The muon-neutrino and muon-antineutrino mean free paths as a function of the ratio between 
nucleon and nuclear saturation densities. Results for neutrinoless matter are shown in upper panels, whereas 
results for neutrino trapping with Yi = 0.3 are shown in lower panels. (Right) Same as in the left panel, but for 
electron-neutrino and electron-antineutrino [34]. 



neutron star. If this happened, we could say that the neutrino is trapped inside the neutron star. The 
"neutrino trapping" can have a serious consequence in the properties of neutron stars. Indeed, it plays an 
important role in the supernova explosion or the neutron star cooling process. 

Although in the Standard Model massless neutrinos have zero magnetic moment and electronic 
charge, there are some evidences that neutrino-electron has a magnetic moment (p v ) smaller than 1.0 x 
lO -10 /^ at 90% confidence level [30], where n# is the Bohr magneton. A stronger bound of /u v < 
3.0 x lO -12 ^ also exists from astrophysical observation, particularly from the study of the red giant pop- 
ulation in globular clusters [31]. On the other hand the neutrino charge radius of v e has been estimated 
from the LAMF experiment [32] to be R 2 = (0.9 ± 2.7) x 10~ 32 cm 2 = (22.5 ± 67.5) x 10" 12 MeV~ 2 , 
while the plasmon decay in the globular cluster star predicts the limit of e v < 2 x 10~ 14 e [31], where e 
is electron charge. 

The NMFP (symbolized by X) is obtained by integrating the cross section of the neutrino-dense- 
matter interaction over the time- and vector-component of the neutrino momentum transfer, 

1 f^Ey-qo r2E v 2 %\q\ 1 (Pa 



X(E V ) J qo ni J ^ u E^E V V d 2 D.'dEy 
where the cross section for every type of neutrino can be calculated from 



1 d 3 o G F E' 



v Im(L„ v IF v ). (38) 



V d 2 CUdE' v 32n 2 E v 

where L piV is the neutrino tensor and IT^ V is the target polarization tensor, which defines the particle 
species. Details of the procedure for calculating Eqs. (I37T ) and (I38T ) can be found in Ref. [33] 

In order to see the effect more clearly, in this calculation we use relatively large values of neutrino 
dipole moments and charge radii. Furthermore, besides using conventional matter, we also consider mat- 
ter with neutrino trapping. The existence of the neutrino in matter allows for the presence of a relatively 
large number of protons and electrons compared to the case of neutrinoless matter. The appearance of 
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Fig. 14. Total differential cross sections as functions of the momentum transfer (f = (qo,qi) for the case where 
neutrinos are trapped with F/ £ = 0.3 and baryon density p# = 5po, while E v = 5 MeV [35]. 

these constituents is then followed by the appearance of a small number of muons at a density larger than 
two times the nuclear saturation density. 

The results for muon-neutrino and electron-neutrino mean free paths are shown in Fig. [T3] where we 
chose the neutrino energy E v = 30 MeV. From this figure it is obvious that in the case of neutrinoless 
matter, contribution from the neutrino charge radius is too small compared to that from the nucleon weak 
magnetism. On the other hand, in the case of neutrino trapping, but with zero muon-neutrino dipole 
moment and charge radius, the difference in the neutrino and antineutrino mean free paths is significant, 
especially at high densities. Although almost similar to the neutrinoless matter case, contribution from 
the neutrino charge radius is very small, and therefore it is invisible in this kinematics. It is also apparent 
that in both cases for densities around (2 — 3)po, and X v behave differently, i.e., if neutrinos are present 
in matter, the mean free paths increase the density, but if neutrinos are absent, the opposite phenomenon 
is observed. We have shown that this result is more pronounced at higher neutrino energies [34]. 

In the case of electron-neutrino, although the effect is less significant for neutrinoless matter at high 
densities, contribution from the neutrino charge radius yields an enhancement to the difference between 
X y and Xy. On the contrary, for the case of zero neutrino dipole moment and charge radius, but with 
neutrino trapping, the difference of neutrino and antineutrino mean free paths is suppressed. Thus for 
this kind of matter, we may conclude that the neutrino charge radius does not contribute to the mean free 
path difference. 
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Fig. 15. Critical densities of different RMF models as a function of the neutrino fraction in matter. 



Another alternative to study the properties of the neutrino-dense matter interaction is by investi- 
gating the sensitivity of the corresponding differential cross section to these possibly nonzero neutrino 
electromagnetic properties. The result for the case of neutrino trapping with lep tonic fraction Y\ e = 0.3 
and baryon density p« = 5po is shown in Fig.[l4j where we plot the total differential cross sections as 
functions of the time- and vector-component of the four-momentum transfer cf. 

From Fig. [14] we can clearly see that the cross section magnitude is quite sensitive to the neutrino 
electromagnetic properties. Including contribution from the neutrino magnetic moment and charge radius 
leads to a cross section which is 14 times larger than that of point neutrino. It is also apparent that most of 
this contribution come from the neutrino charge radius, whereas the different shape in the q\ distribution 
is mainly driven by the neutrino magnetic moment. 

We have also compared this result to the case of neutrinoless matter [35]. We found that, in general, 
neutrino trapping enhances the magnitude of the cross section, as we expected, since matter with neutrino 
trapping contains more protons and, hence, electromagnetic interactions are more likely to occur. If we 
switched off the neutrino electromagnetic properties, this behavior clearly disappears. 



3.4. Critical Density of Dense Matter 

From Fig. [TT] it is apparent that a neutron star is expected to have a solid inner crust above its liquid 
mantle. The mass of this crust has been known to be sensitively dependent on the density of the inner 
edge and on the EOS. Therefore, investigation of the dynamical instability of dense matter will certainly 
shed important information on the inner structure of neutron stars. 

Recent investigations on this topic have disclosed that the critical density (p c ), a density at which the 
uniform liquid becomes unstable to a small density fluctuation, can be used as a good approximation of 
the edge density of the crust [36]. The important aspects of their finding are that the inner edge of the 
crust density to be p e d ge = 0.08 fm~ 3 and a measurement of the neutron radius in the 208 Pb nucleus will 
provide useful information on the p c . 

To this end we have investigated the behavior of the predicted critical density (p c ) in matter with 
and without neutrino trapping for several different relativistic mean field (RMF) models [37]. The ratio 
between the critical density and the normal nuclear density as a function of the neutrino fraction in matter 
is shown in Fig.[l5]for several different RMF models. 

It is found that different treatments in the isovector-vector sector of RMF models yield more sub- 
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Fig. 16. Connecting and extending from what we have learned so far. 



stantial effects in matter without neutrino trapping than in matter with neutrino trapping. Moreover, for 
matter with neutrino trapping, the value of p c does not significantly change with the variation of the 
models nor with the variation of the neutrino fraction in matter. In this case, the value of p c is larger for 
matter with neutrino trapping. These effects are due to the interplay between the major role of matter 
composition and the role of the effective masses of mesons and nucleons. 

It is also found that the additional nonlinear terms of Horowitz-Piekarewicz and Effective RMF 
models can prevent another instability at relatively high densities from appearing. This can be traced 
back to the effective a mass, which goes to zero when the density approaches 6.5po- 

4. Connecting What We Have Learned So Far 

So far, we have discussed hypernucleus and neutron star starting from the Pauli exclusion principle. 
We know that the important ingredient of a hypernucleus is hyperon, whereas the basic building block 
of a neutron star is neutron. According to their symmetry, both neutron and hyperon are classified in 
the SU(3) baryon octet. We have also implicitly built an argument that without this principle we know 
nothing about the physics of hypernuclear and neutron star. 

Now, we can nicely connect hyperon, neutron, hypernucleus, and neutron star as shown in Fig. [16] 
Starting from this figure, we may ask: what can be constructed based on our knowledge about them? A 
tentative answer is also shown in this figure. 

Within the RMF model, the question whether a neutron star could be considered as a giant hypernu- 
cleus has been asked more than 20 years ago [38]. Earlier discussion which showed that the existence 
of hyperon in neutron star is plausible has been done based on Fermi gas model [39]. For our current 
discussion this would not be too difficult to understand if we look back at Fig. (TTJ where we can clearly 
see that at certain density (depth) there is a A population. In fact, Glendenning has found that this 
population starts to grow at baryon density p# ~ 0.4 fm _1 and starts to dominate all star ingredients at 
Pb ~ 0.9 fm _1 [38]. Such population is found in the cores of the relatively heavier neutron stars and is 
about 15%-20% of the total baryon population. 

We can continue the discussion to the hyperon stars as well as quark stars. They are indeed more 
massive, more complicated, and more interesting than the neutron star. However, we have to stop this 
discussion here, otherwise I have the feeling that our experimentalist colleagues would complain that we, 
theoreticians, are too imaginative and too creative. Could this be true? The next section could convince 
you that this could be true. 
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5. Econophysics, are Theoreticians too Creative? 

5. 1 . Why Econophysics ? 

The scientific challenge to understand the nature of complex systems is too tempting to the physicists. 
Some of them have been fascinated by the huge, and also growing, amount of economics data recorded 
minutes by minutes for decades. Among these interesting data, the fluctuation of stock exchange indices 
is of special interest, since it might indirectly reflect the economic situation in a certain region. Fur- 
thermore, the advancement in computing capabilities has enabled us to handle this large amount of data, 
unlike the situation almost 40 years ago, when Mandelbrot investigated approximately 2000 data points 
of cotton prices [40]. 

It is said that such studies could explain the nature of interacting elements in the complex systems 
and, therefore, could help to forecast economic fluctuations in the future. In other words, these studies 
were intended to produce new results in economics, which might help us to avoid economic "earth- 
quakes" such as what happened in Indonesia several years ago [41]. 



5.2. Time Series Analysis 
We start our analysis with an investigation of the daily index returns, which are defined as 

Z At (t) = lnY(t + At)-lnY(t) = \n[Y(t + At)/Y(t)], (39) 

where Y(t) indicates the closing index of the stock at day t. We use the time series data from two different 
stock indices, the Jakarta Stock Exchange Index (abbreviated with IHSG, an acronym of Indeks Harga 
Saham Gabungan or the composite stock exchange price index) and the Kuala Lumpur Stock Exchange 
index (KLSE), which belong to different countries. Comparing the two indices would be very interesting 
since both Indonesia and Malaysia underwent the same monetary crisis in 1997, which was then followed 
by financial crashes in almost all economic sectors, but with quite different economic situations. 

Figure [FT] shows the time series of the IHSG and KLSE indices along with their logarithmic returns 
calculated by using Eq. (|39l >. For comparison, in this figure (top panels) we also display the historical 
time series of the exchange rate of Indonesian Rupiah and Malaysian Ringgit, since the economic crashes 
stalled with the devastating decline of these rates. It is obvious from this figure that the fluctuation Y(t) 
of both indices is more dramatic than that of S&P 500, indicating that in this case the situation is far 
more complex. After the crash the magnitude of returns Z(t) is significantly larger in both indices, or, in 
the economics language, the probability to gain or to loose becomes larger than before. 



5.3. Omori Law, from Geophysics to Finance Physics 

To accurately determine the time position of the crash we make use of the Omori law. Omori law, which 
was originally used in geophysics [45], states that the number of aftershock earthquakes per unit time 
measured at time t after the main earthquake decays as a power law. Practically it is written as 

n{t) = K(t + x)-P, (40) 

where K and x are two positive constants. The cumulative number of aftershocks is obtained by integrat- 
ing Eq. (|40l from to t, i.e., 

w = /Vgy-n-' - *["+<>'-"-^] . (4i, 

Jo 1 - p 

This number is found to be more relevant to compare with the presently analyzed data. By calculating the 
number of returns given by Eq. (l39l that exceeds one standard deviation (a) up to 7a, we arrive at Fig.fTFI 
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Fig. 17. (Left) The ratio between 1000 Indonesian Rupiah and 1 US Dollar (top), the Jakarta Stock Exchange Index 
(IHSG) Y(t) (middle) and the logarithmic return Z(t) calculated from Eq. (|39l as a function of time t sampled 
with Af = 1 day (bottom). (Right) The ratio between 1 Malaysian Ringgit and 1 US Dollar (top), the Kuala 
Lumpur Stock Exchange Index (KLSE) Y(t) (middle) and the logarithmic return Z(f ) (bottom) as a function of 
time sampled for Af = 1 day. On 1st September 1998 Malaysia imposed currency controls, the Ringgit was pegged 
with US Dollar with a fixed rate of 3.80 Ringgit per Dollar. The arrows in the bottom panels indicate the time 
position when the IHSG and KLSE started to crash. Data are taken from Ref. [42] (IDR/USD and RM/USD 
exchange rates) Ref. [43] (IHSG stock index), and Ref. [44] (KLSE stock index). 




t (day) 

Fig. 18. Number of the returns given by Eq. d39l that exceeds one up to seven standard deviations (a) for the IHSG 
(top panel) and the corresponding derivatives at la (lower panel). 
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Fig. 19. (Left) Log-log plots of the IHSG probability of return to the origin as a function of the sampling time 
Af before and after crash. Solid lines are obtained from linear regressions to this probability. The slopes of these 
lines, which are equal to negative inverse of the Levy stable distribution indices, are shown in the left panels. 
Comparison between normalized distribution functions for Af = 1 and normalized Levy stable distributions with 
parameters obtained from the left figures is shown in the second column panels. The standard deviation and the 
number of data used are also shown in these panels. (Right) Same as in the four left panels, but for the KLSE case 
before and after crash. 



for the IHSG case. From this figure we can clearly see that there exists a discontinuity on 5 August 
1997, which indicates the start of a crash process. To exactly find the position of this discontinuity we 
calculate the first derivative of N(t) for N(t) obtained with la, shown in the lower panels of Fig. [TU 
From these three panels we are convinced that the IHSG started its crash process on 5 August 1997. The 
same procedure has been also applied to the KLSE index, from which we obtain that the KLSE started to 
crash on 10 July 1997, almost one month before the IHSG case. The crash positions of the two indices 
are indicated by the vertical lines in Fig.fTTl 



5.4. Probability Density Function 

Following previous studies [46,47] we investigate the probability density function [PDF or P(Z)] of the 
return to the origin P(0) in order to investigate the statistical properties of both indices. The advantage of 
using such analysis is that we can reduce the statistical inaccuracies, since the number of data included 
is relatively small, whereas the probability is largest at Z = 0. 
Starting with the characteristic function [47] 

q>(Af) = e-^ a , (42) 

the Levy stable distribution is given by 

P(Z) = - e-^ a cos(qZ)dq . (43) 



71 JO 

From Eq. d43l the probability of return to the origin P(0) reads 



P(0) = QVg) , ( 44) 
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Fig. 20. Moments of the distribution of the IHSG (top) and KLSE (bottom) normalized returns, before (left panel) 
and after (right panel) crash, given in Eq. d45b for At = 1,2,3,4,6,9, 13, and 19 days compared with those obtained 
from the Gaussian distribution (solid thick lines). 



where F indicates the Gamma function. 

The log-log plots of P (0) as a function of the sampling time At for both IHSG and KLSE indices 
are shown in the first and third column panels of Fig. [19] The slopes of linear regressions to these plots 
equal the negative inverse of the Levy stable distribution indices a. Using this index we calculate the 
parameter yby means of Eq. (l44l) and plot the "theoretical" PDF as a function of normalized returns Z/a, 
using Eq. (l43l) . where a is the standard deviation of the distribution, and compare it with the empirical 
PDF obtained from data in Fig. [19] Figure[l9] reveals the fact that in the case of IHSG the distribution 
turns from Levy to Gaussian after the crash. Surprisingly, the KLSE stock index shows a contrary result, 
the distribution alters from a Gaussian to a Levy one after the crash. 



5.5. Distribution Moment 

It has been argued that the use of the return probability to the origin PDF(0) to estimate the Levy stable 
distribution index a is statistically not optimal, due to discreteness of the distribution [48]. Instead of 
exploiting such method, a different strategy has been proposed, i.e. by calculating a by means of the 
slope of the cumulative distribution tails in a log-log plot. To further test their results on the scaling 
behavior, Refs. [48,49] analyzed the moments of the distribution of normalized returns 

m = (\ g (t)\ k ). (45) 

where the normalized returns g(t) is defined by 

z(0-(z(0) 



V(z 2 (0) - (z(0) 2 



(46) 



with (Z(t)) the time average of Z(f) over the entire of time series. In the case of the S&P 500 index 
the result is found to be consistent with the analysis of the tails of cumulative distributions. References 
[48,49] pointed out that the change in the moments behavior originates from the gradual disappearance 
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of the Levy slope in the distribution tails. In our case it is also important to cross-check the result shown 
in Fig.[l9j since the number of data used in this analysis is significantly smaller than that of previous 
analysis on the S&P 500 index, e.g. Ref. [48], which used approximately 5 x 10 6 data points. 

It has been also shown in Ref. [48] that Eq. (1431 ) will diverge for k > 3. In this study we also constrain 
k within < k < 3. The results for both indices compared with the moment obtained from a Gaussian 
distribution have been reported in Ref. [50]. It is found that the IHSG moment retains its scaling up to 
At = 9 days, only after At = 13 the moment starts to deviate toward the Gaussian distribution. In the 
KLSE case the moment quickly converges to the Gaussian distribution as soon as At increases from 1 
day and does not show any scaling behavior as in the former case. 

In Fig.[20] we display the moments of both stock indices in the case of before and after the financial 
crash. A consistent result, compared to the analysis of Fig.[[9l is obtained from Fig.|20l i.e., in the case of 
IHSG the distribution becomes closer to Gaussian after crash, whereas the KLSE moments move away 
from the Gaussian distribution after the crash. 

Another important finding obtained from Fig. [20] is that the scaling behavior up to At = 9 shown by 
the IHSG case in the previous section originates from the "before crash" period. After the crash, the 
IHSG moments quickly converge to a Gaussian distribution. In fact, this phenomenon has already been 
seen in the left panels of Fig.QjS where the empirical probability of return to the origin is more scattered 
after crash. 
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